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SYSTEM AND METHOD FOR PERIODICALLY 
ADAPTIVE GUIDANCE AND CONTROL 

BACKGROUND OF THE INVENTION 

Field Of The Invention 

[0001] The present invention relates to guidance and control systems, and more 
particularly, to adaptive guidance and control systems. 

Background Information 

[0002] When guiding a device, such as an interceptor missile, toward a 
moveable object or target, it is known to address target maneuvering in the 
control of the missile. For example, the object may be configured to induce a 
miss distance by initiating an acceleration maneuver. Against a missile utilizing a 
proportional guidance system, target maneuvering can be controlled as a function 
of any known missile parameters, such as the missile's autopilot time constant (its 
agility) and effective navigation ratio. If the target does not have a priori 
knowledge of the missile parameters, does not know time to intercept, and/or does 
not have visual or electronic contact with the missile, then periodic maneuver 
sequences such as a barrel roll or weave maneuver can be executed. In addition to 
calculated maneuvering, targets such as tactical ballistic missiles (TBM) can 
unintentionally spiral or weave into resonance as they re-enter the atmosphere due 
to either mass or configuration asymmetries. 



[0003] The control of a moving object to a maneuverable target is discussed in 
the book, "Tactical and Strategic Missile Guidance", Third Edition, Volume 176, 
by Zarchan, and in a Zarchan article, "Proportional Navigation and Weaving 
targets", Journal of Guidance, Control, and Dynamics, Volume 18, No. 5, 
September-October 1995, pages 969-974, both of which are incorporated herein 
by reference. Zarchan discloses a scalar classical pro nav guidance law as: 
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wherein the missile acceleration command (n c ) is perpendicular to the target 
line-of-sight vector for small line-of-sight angles (X ) as evidenced by the X = (y 
/ r^) approximation. Small angle approximations are used to simplify modern 
guidance laws that use scalar Y axis Cartesian Kalman filter estimates of 
target-to-missile range and velocity (y and y). For a dual-range air-to-air missile, 
however, guidance laws with these small line-of-sight angle approximations result 
in excessive missile divert requirements, which waste valuable on-board energy, 
and can degrade performance. 

[0004] "Line of Sight Reconstruction for Faster Homing Guidance", presented 
as Paper 83-2170 at the AIAA Guidance and Control Conference, Gatlinburg, 
Tenn., August 15-17, 1983; received August 24, 1983: revision received January 



3, 1984. Copyright ® 1983, by F.W. Nesline and P. Zarchan, which is 
incorporated herein by reference, describes line-of- sight ( LOS) reconstruction. 
This Zarchan document discloses another scalar pro nav guidance law as: 

n = ^ )+ZEM ™ } - C A n L (2) 

where ZEM m = 0.5^ 2 o n T (3) 

[0005] The missile acceleration command (n c ) is over amplified by N' and 
negated by an achieved missile acceleration quantity (C 4 nj to make-up for 
expected dynamics lag in the missile's autopilot. Zarchan also describes a 
weaving target guidance law with a priori estimates of the target weaving 
frequency. Equation (2) uses an estimated target acceleration (n T ) presented in 
Equation (3) as the zero effort miss (ZEM^) due to a target maneuver. 
[0006] Zero effort miss is defined as the distance the missile would miss if the 
target continued along its present trajectory and there were no more missile 
acceleration commands. A condition of Equation (3) is that the zero effort miss is 
caused by constant target acceleration for all time-to-go before intercept. Zarchan 
derives the zero effort miss due to a target maneuver frequency (co T ) with target 



acceleration and acceleration rate magnitudes (n x and ri T ), and documents the 
scalar weaving target guidance law as: 
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[0007] However, known systems do not address how to estimate a target 
maneuver frequency (oo T ) for the weaving target guidance law, and use small angle 
approximations. 

[0008] U.S. Patent No. 4,494,202 (Yueh), the disclosure of which is hereby 
incorporated by reference, describes a fourth order predictive augmented 
proportional navigation system terminal guidance design with missile attached 
target decoupling. In Yueh's system, the estimated value of target lateral 
displacement is derived from a representation of the measured LOS angle 
combined with radome error and noise in the system, and small angle 
approximations are used. The gain C 3 associated with a predetermined estimate of 
target maneuver augments the target acceleration term into a proportional 
navigation design. However, this gain is not adaptive, and relies on the target 
bandwidth (i.e., a predetermined value based on expected target maneuver). 
Thus, Yueh discloses an angular guidance approach which allows guidance for all 



line-of-sight (LOS) angles, yet maintains the conventional use of small angle 
approximation. Yuen does not support the framework for building a fully 
descriptive target state equation for adaptively estimating a target maneuver 
frequency in Cartesian space. 

[0009] Warren, Price, Gelb and Vander Velder in "Direct Statistical Evaluation 
of Nonlinear guidance systems" AIAA Paper No. 73-836, which is incorporated 
herein by reference, disclose several missile guidance laws. In particular two 
optimal linear guidance laws C and D are shown that use Kalman filters which are 
synonymous with the classical optimal guidance control laws of Zarchan. 
Guidance law C minimizes the performance index, neglecting autopilot dynamics, 
resulting in C 4 =0 and C 3 is a function of the target bandwidth, A t . Guidance law 
D is derived including the first-order autopilot model dynamics and the resulting 
n' and C 4 are both functions of A m and t g(r The gains and navigation ratios for the 
Classical optimal guidance control laws of Warren et al., and Yueh are shown in 
Table 1. 
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where: X T =A T t go , X M = X M t go ■ t go = Estimated Time-to-Go [sec] 

X T = Target Maneuver X M = Missile Autopilot/ Airframe y = Optimal quadratic Cost 

Bandwidth [rad/sec] Bandwidth [rad/sec] Function Weighting Factor 



[0010] As with the systems disclosed by Zarchan, Yueh's system and those 
disclosed in Warren, are subject to the classical augmented proportional navigation 
problem, the zero or constant maneuver frequency assumptions, thus wasting 
valuable energy against highly maneuverable targets. 

SUMMARY OF THE INVENTION 
[0011] The present invention is directed to a system and method for adaptively 
guiding a moving device toward a maneuverable object. Exemplary embodiments 
include means for generating a guidance command signal from a vectored line-of- 
sight (LOS) between a device and an object using a position parameter of the 
object relative to a guidance frame, and an estimated object state produced in the 
guidance frame using the vectored line-of-sight; and means for transmitting the 
guidance command signal to an on-board guidance control of the device. 
[0012] Exemplary embodiments include means for generating a signal 
representing a predicted position of the object from object position parameters 
relative to a guidance frame and periodically adaptive estimated object state 
produced in the guidance frame using the object position parameters; and, means 
for transmitting the signal to an on-board guidance control of the device. 



BRIEF DESCRIPTION OF THE DRAWINGS 
[0013] These and other features and advantages of the claimed subject matter 
will become more apparent from the detailed description and accompanying 
drawings, in which: 

[0014] Figure 1 is a periodically adaptive all angle guidance system in 
accordance with an exemplary embodiment of the present invention; 
[0015] Figure 2 is representation of normalized target acceleration power 
spectral density of a harmonically balanced Kalman filter bank in accordance with 
an exemplary embodiment of the present invention; 

[0016] Figure 3 is a block diagram of Vectored LOS reconstruction for a 
Gimbaled seeker in accordance with an exemplary embodiment of the present 
invention; 

[0017] Figure 4 is a block diagram of Vectored LOS reconstruction for an 
electronically steered seeker in accordance with an exemplary embodiment of the 
present invention; 

[0018] Figure 5 is a block diagram of a periodically adaptive guidance system 
in accordance with an exemplary embodiment of the present invention; 
[0019] Figure 6 is a block diagram of a three dimensional periodically adaptive 
target state estimation system; 

[0020] Figure 7 is a block diagram of a harmonically balanced Kalman filter 
bank in accordance with an exemplary embodiment of the present invention; 



[0021] Figures 8a-8b are representative paths and power spectral density charts 
for a missile intercepting a target according to Classical Optimum Control and 
Periodically Adaptive Control, respectively, in accordance with exemplary 
embodiments of the present invention; 

[0022] Figures 9a-9d are exemplary response capability diagrams produced 
using periodically adaptive all angle guidance for a periodic maneuver target, a 
constant acceleration target, a time varying jinking target, and a constant velocity 
target, respectively; 

[0023] Figure 10 is a block diagram of a guidance system with LOS and range 
vectors in accordance with an exemplary embodiment of the present invention; 
[0024] Figure 1 1 is a representation of a fire control guidance system for a 
missile in accordance with an exemplary embodiment of the present invention; and 
[0025] Figure 12 is a block diagram of periodically adaptive guidance for fire 
control in accordance with an exemplary embodiment of the present invention. 

DETAILED DESCRIPTION OF INVENTION 
[0026] Figure 1 shows a guidance system 400 for guiding a device toward an 
object according to an exemplary embodiment of the present invention. 
[0027] The system 400 includes means 460 for generating a guidance command 
signal from a vectored line-of-sight (LOS) between a device and an object using a 
position parameter of the object relative to a guidance frame, and an estimated 
object state produced in the guidance frame using the vectored line-of-sight, and 
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means 463 for transmitting the guidance command signal to an on-board guidance 
control of the device. The guidance command signal generating means can be 
configured using a single processor, computer or other electronic means, or via 
any number of processors, computers, etc., or combination thereof for 
implementing the function. The guidance command signal generating means can 
include means for creating a vectored line-of-sight (LOS) between the device and 
the object using the position parameter, generally represented as 461 . The means 
for creating a vectored LOS can include LOS reconstruction. Line-of-sight 
reconstruction (LOS) provides an ability to have a fast missile autopilot time 
constant without causing guidance loop destabilization because of the relatively 
lower bandwidth of the seeker tracking control loop. Vectored LOS 
reconstruction is direction cosine rate integration of the seeker gyro rates and 
vector addition of the tracking errors. The result is a reconstructed target 
line-of-sight vector that is accurate for all LOS angles, thus creating an object line- 
of sight (LOS) in a guidance frame. Vectored LOS reconstruction combined with 
missile (IMU) information 450 transforms the target range, seeker gyro rates and 
tracking error measurements into a target-to- missile range measurement vector 
with components in a three-dimensional inertial guidance frame. 
[0028] Estimated object states are produced relative to the guidance frame 
using a sequential object line-of-sight LOS. The estimated object state producing 
means is a periodically adaptive object state estimator represented by block 464, 
can include object state estimators with Kalman filters. Three-dimensional object 



state estimators are decoupled into three single dimension object state estimators 
that each operate independently on one of the range measurement vector 
components. The linear independent models describe the target-to-missile 
dynamics in each guidance frame axis. The three uncoupled object state 
estimators quickly and accurately adapt to the object acceleration frequency 
observed in each guidance frame axis by using a harmonically balanced multiple 
model Kalman filtering algorithm, as disclosed in "Multiple Model Adaptive 
Estimation and Prediction With The Harmonically Balanced Kalman Filter Bank", 
A Thesis Presented To The Faculty Of California State Polytechnic University, 
Pomona, by Randall Joseph Muravez, 1989 which is incorporated herein by 
reference. In addition, the three uncoupled periodically adaptive target state 
estimators and the multiple model Kalman filters within them can be executed in 
parallel to maximize processor throughput (i.e. digital array processing). 
[0029] In the system 400, an interdependent group of variables (states) work 
together as a unified whole. The states of a system can be arbitrarily defined. 
The object states utilized in exemplary embodiments of periodically adaptive all 
angle guidance include range, velocity, acceleration and a further time derivative 
(i.e. acceleration rate). 

[0030] The Means for creating guidance command signals using the estimated 
object state can include the implementation of various guidance and control laws. 
All Angle Guidance and Control 465 can include guidance command laws 
including All Angle Classical Pro Nav Guidance, Periodically Adaptive Guidance 
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and All Angle Augmented Pro Nav Guidance which are applied to create guidance 
command signals perpendicular to the target LOS vector for all LOS angles, in 
order to null the angular rate of the target LOS vector with the missile acceleration 
commands. 

[0031] A missile or device acceleration guidance command signal is fed to the 
missile autopilot 456 via wire or wireless means 463. The missile autopilot 456 
communicates with the missiles inertial measurement unit (IMU) 450 as well as a 
propulsion system and/or control mechanism (i.e., aerodynamic control surfaces) 
of the interceptor (i.e., missile, device). 

[0032] Figure 1 also shows means for obtaining an object position parameter, 
generally 425. In an exemplary embodiment object position parameters are 
determined from target to seeker range and angular measurements. 
[0033] "All Angle Guidance" means that for all possible target-to-missile 
line-of-sight angles, the guidance algorithms will perform accurately because 
vector algorithms are used to eliminate small angle approximations. Applying 
"All Angle Guidance" vector algorithms to current algorithms leads to higher 
probability of kill and lower missile agility requirements for all types of 
engagement geometries. An example of where "All Angle Guidance" vector 
algorithms can be required is for a short-range "dogfight 'engagement where, the 
close proximity of the object and the interceptor create a dynamic environment 
were small angle approximations can be of no value. These algorithms can also be 
beneficial for the final aim-point portion of longer range engagements. 
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[0034] An all angle classical proportional navigation (pro nav) guidance law is 
derived by applying "All Angle" vector algorithms to Zarchan's classical pro nav 
guidance law of Equation (1) so that the missile acceleration command vector (nj 
is always perpendicular to the target line-of-sight unit vector (I T ): 
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[0035] The closing velocity (V c ) and approximated time-to-go (t go ) are defined 



as: 
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and the "All Angle" line-of- sight rate vector (X) to be nulled is defined as: 
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[0036] The navigation gain (N') for the classical pro nav guidance law is 
typically a constant. The target-to-missile range and velocity vectors ( and 
Ytm) come from combining estimates from the three uncoupled target state 
estimators. The line-of-sight rate vector (XX of the "All Angle" classical pro nav 
guidance law near the left-hand side of Equation (5) illustrates more intuitively 
that {il is to be nulled, and does not use the time-to-go approximation. 

An all angle optimal augmented pro nav guidance law is derived from 
applying "All Angle" vector algorithms to Zarchan's scalar modern optimal 
augmented pro nav guidance law of Equation (2) so that the missile acceleration 
command vector (nj is always perpendicular to the target LOS unit vector I x : 
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[0037] The relationship shown in Equation (4) is used to keep the classical pro 
nav part of the new law intact as function of X , which is more intuitive and does 
not use the time-to-go (t go ) approximation of Equation (3). Zarchan's modern 
navigation gain (N') and the achieved missile acceleration quantity (C 4 nj that 
optimally make-up for the expected dynamics lag in the missile autopilot remain 
unchanged. The classical guidance augmentation using estimated target 
acceleration (n T ) is presented in Equation (11) as the zero effort miss vector 
(ZEMtm) caused by a constant target maneuver. An important enhancement of the 
new law in Equation (10) is the application of two cross-products to ZEMtm using 
the target line-of-sight unit vector (I T ). This results in properly applying the 
ZEMtm augmentation so as to always be perpendicular to the target LOS unit 
vector I x . 

[0038] Thus, the "All Angle Guidance" law of Equation (10) can lead to higher 
probability of kill and lower missile agility requirements for all types of 
engagement geometries. 

[0039] Another feature used in accordance with exemplary embodiments is 
"Periodically Adaptive Guidance." This means that the target acceleration 
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frequency (spiraling, weaving, none) is periodically adaptively estimated and used 
to successfully improve missile guidance. The result can lead to higher probability 
of kill and lower missile divert requirements for a wide range of maneuvering 
targets. 

[0040] A periodically adaptive guidance law is derived from improved 
prediction of ZEM m based on probability weighted multiple models of different 
target maneuver frequencies and their target acceleration and acceleration rate 
states as described for one dimension in "Multiple Adaptive Estimation and 
prediction with the harmonically balanced Kalman Filter Bank" : 
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where: 



0 M = atan2(o> Tm ,X Tm ) 



= X Tm +u Tm 



®Bm = d-(4^ / C0jj) / COL 
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[0041] As indicated in Figure 1, the multiple model probability weighting 
vectors (WJ defined inequations (12) and (13) originate from the periodically 
adaptive target state estimators that are based on multiple model Kalman filtering. 
The probability weighting vectors (W^ indirectly represent a three-dimensional 
estimate of the target maneuver frequency oo t and its correlation time constant T Cm . 
The probability weighting vectors (WJ contain weights , wherein the sum of the 
probability weights for any axis of the guidance frame is unity, and they uniquely 
make the "Periodically Adaptive Guidance" law of Equation (12) periodically 
adaptive. 

[0042] As shown in Equation (12) using the vector definitions of Equations 
(13) and (14), the m* Kalman filter model's estimated target acceleration and 
acceleration rate state vectors (A T m and Arm) are multiplied by the m* model's 
state transition elements $ 13ra (t go ) and <£ 14ra (t go ). Then they are multiplied by the 
probability weighting vectors W m and summed over the total number of models in 
the Kalman filter bank (NMB). The m" 1 model's state transition elements 0 13m (t g0 ) 
and 3> 14m (t go ) are time-to-go predictive coefficients based on an m 111 model's a priori 
hypothesis of expected target maneuver frequency ( co^) and its correlation time 
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constant (T C J. The state transition elements are presented in equations (13) and 

(14) , their derivation is shown in "Multiple Model Adaptive Estimation and 
Prediction with the harmonically Balanced Bank". 

[0043] The m Ul model state transition elements are obtained using equations 

(15) through (17), with the variables A Xm and o> Tm established through 
parameterization of the harmonically balanced Kalman filter bank. To periodically 
adaptively estimate states of a maneuvering target, it is desired to adapt to a 
constant acceleration target (T c * T c max [Sec], co M = 0 [R/S]), to a constant 
velocity target (T c * 0.2 [Sec], o) M = 0 [R/S]), and to a jinking or other type of 
maneuvering target (T c = unknown [Sec], o) M = unknown [R/S]). The common 
parameterization design approach is to determine each Kalman filter's dynamical 
parameters by discretizing a continuous parameter space, which is accomplished 

by dispersing a set of dynamical parameter vectors [d x d^g] throughout a 

region of reasonable parameter vector values. 

[0044] Figure 2 shows an exemplary harmonically balanced Kalman filter bank 
model 200. The parameterization of the "harmonically balanced Kalman filter 
bank," shown in Figure 1, is accomplished with the m 01 periodic dynamical 
parameter vector d m shown by Equation 18: 
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where A Tm = m ai model target bandwidth, and u) Tm = m ai model target maneuver 
frequency. 

[0045] A hypothesis set of periodic dynamical parameter vectors is split into two 
parts: 

KiW A Wi d NMB~l nQ 

where 

[0046] The first dynamical parameter vector d! defines the highly time- 
correlated "zero-frequency" system dynamics hypothesis model. For Periodically 
Adaptive Guidance Control, d m at m=l defines the constant target acceleration 
hypothesis, shown in equation (20): 
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where T C MAX [Sec] is the maximum periodic system dynamics time constant. 
The remaining set of the periodic parameter vectors [d J 2:NMB encompass a periodic 
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parameter space 201 bounded by T C MAX, shown as line 202 and the maximum 
system dynamics maneuver frequency m T MAX represented as line 203. The 
periodic parameter space 201 with these maximum boundaries is the domain of the 
normalized PCA power spectral density function. Each periodic parameter vector 
included in [dJ 2:NMB has an identical balanced maximum time constant T C MAX 
[Sec] or damping coefficient, which is equal to 2/T c MAX [1/Sec]. And, each 
periodic parameter vector is harmonically separated by a constant half-power 
bandwidth: \!T C MAX [Rad/Sec]. Therefore, the parameterization design of the 
"harmonically balanced Kalman filter bank" is completed with the following set of 
periodic dynamical parameter vectors: 
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[0047] The total minimum number of parameter hypothesis models in the 
"harmonically balanced Kalman filter bank" (NMB) is simply a "periodic 
parameter space area" bounded by the maximum periodic dynamical parameters 
T C MAX [Sec], and <s> T MAX [Rad/Sec], plus one hypothesis model shown as line 
204 based on the "zero-frequency" dynamical parameter vector d,,: 
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NMB = <s> T MAX T C MAX + 1 
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and models based on d m ;m=2 to NMB are represented generally in Figure 2 as 
waveforms 206. 

[0048] The maximum periodic parameters T C MAX and <x> T MAX are specified 
based on a designer's discretion for a particular application. As a result, many 
arrangements of the "harmonically balanced Kalman filter bank" are possible from 
the laws of Eqs. (20) and (21). 

[0049] The minimum time constant T C MIN, shown as line 205, is defined as 
the height of an "equally weighted" harmonically balanced PCA spectral density 
function, which is a band limited white noise acceleration or constant velocity 
model with minimum correlation time constant T c MIN, similar to an ideal low- 
pass filter. The minimum system dynamics time constant for the "harmonically 
balanced Kalman filter bank" is inversely proportional to the maximum system 
dynamics maneuver frequency: 



T C MIN * 



{xi T MAX 



[Sec] 



(23) 
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[0050] The other maximum periodic parameter, the maximum maneuver 
frequency id T MAX, is bounded by the maximum realistically expected target 
maneuver frequency for a particular target based on a priori knowledge. In 
addition, o) T MAX should be maximized according to Eq. (23) so as to minimize 
the target acceleration time constant T C MIN for adequately describing a constant 
velocity target. 

[0051] The "Periodically Adaptive Guidance" law of Equation (12) can 
improve missile guidance. As the estimated target maneuver frequency gets 
higher, the probability weighted state transition elements cause the estimated target 
acceleration and acceleration rate states to have less affect on the zero effort miss 
due to a target maneuver. This is reflected in physics: an acceleration sine wave 
of frequency to that is twice integrated affects its position magnitude by 1/co 2 . 
Thereby, the zero effort miss of Equation (12) used to augment the "All Angle 
Guidance" law of Equation (10) will be smaller for higher target maneuver 
frequencies. For lower target maneuver frequencies the opposite is true, and for 
zero frequency or constant target acceleration, the zero effort miss due to a target 
maneuver becomes the classical augmentation of Equation (11). Additionally for a 
non-maneuvering constant velocity target, the probability weighting vectors can 
become equally weighted and can successfully cancel out the augmentation. 
[0052] Thus, the "Periodically Adaptive Guidance" law of Equation (12) can 
lead to higher probability of kill and lower missile divert requirements for a wide 
range of maneuvering targets. 
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1 . System Integration 

[0053] Figure 1 is a simplified block diagram of a guidance system in 
accordance with an exemplary embodiment of the present invention. Figure 1 also 
shows means for obtaining object position parameters as used in exemplary 
embodiments. 

[0054] The target's actual velocity and acceleration vectors within the guidance 
frame are shown as V G T and A G T and are shown in the target portion of the 
diagram 400. Vectors are generally represented with a capital letter that is 
underlined. The reference frame associated with vectors is a superscript to the 
right of the vector, for example V G is a velocity vector in the G reference frame. 
Measured parameters are indicated with a "*" and subscripts are used to associate 
the parameter with an object, for example V G * is a measured velocity vector 
associated with a target. A over a scalar or a vector indicates an estimated 
value and a " ' " indicates a rate, for example ^ is the estimated change of the 
range with respect to time. 

[0055] The relative geometry of the target in relationship to the missile is 
shown in the relative geometry portion 410 of the diagram where the V G ^ 
represents the closing velocity in vectored form between the missile and the target. 
Subsequently, the range vector are R* 3 ™ represents the end state miss distance 
between the target and the missile in vector form and I G X represents the target LOS 
vector. The target LOS angles k el and k „ are the actual azimuth and elevation 
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angles from the missile to the target. The description of the relative geometry is 
merely a representation of the physical dynamics of the target with respect to the 
guidance frame and the missile. This portion of Figure 1 need not be performed 
by the guidance system, but serve as a representation for illustrative purposes. 
[0056] In the radome portion of the diagram A el and X^ are operated on by 
block 421 [1 + r], where r is the radome slope error. The modified LOS angles are 
combined at 422 with representation of ^ and which are the missile body axis 
angles with respect to the guidance frame, multiplied by r, the radome slope error 
in block 423. The output of the combiner 422 is a measurement of the LOS in 
X\ x and X' ^ which is forwarded to the Gimbaled seeker 430, Radome 
compensation is addressed in U.S. Patent 3,316,549. The Gimbaled seeker track 
loop has Gimbal rates as an input which are the rates of change obtained from a 
rate gyro of the angle between the antenna axis and the missile axis. The gyro 
rates to el *, co^* are integrated in block 432 and the integration is applied to the 
output of the radome A' el and X' a . The combination results in e d , e B , which are 
the angles between the missile to target LOS and the electrical bore sight or track 
errors. The track errors are multiplied by the track loop gain in block 434 and 
that quality along with the gyro rates are applied to the integrator 432 to complete 
the loop. The track errors are also multiplied in block 435 to compensate for the 
radome, resulting in e e i* and e^*. These compensated track errors are inputs to 
the periodically adaptive all angle guidance system. The Gimbal rates in 431 are 
integrated in block 436 and used to compensate the track errors and serve as input 
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to the periodically adaptive all angle guidance system 460. The seeker track errors 
e d * and e^* have radome compensation and target body aim point corrections. 
The gimbal rates, derived from the body rates i|f eZ and are integrated 
resulting in gimbal angles G pitch * and G yaw *, they are applied as inputs to the 
radome compensation 435 and into the periodically adaptive all angle guidance and 
control system 460. An inertial reference unit (IRU) 450, the same as IRU 452, is 
transformed through a transformation matrix C" b into the inertial navigation frame 
and supplied to the periodically adaptive all angle guidance and control system 
460. 

[0057] The output of the periodically adaptive all angle guidance and control 
system 460 results in a missile acceleration guidance command vector in the 
inertial guidance frame n° c±LOS which is transmitted at 463 as an input to the 
missile's autopilot 456. The autopilot generates an acceleration vector which is 
integrated in block 457 to represent a velocity vector. The acceleration and 
velocity vectors are fed into a transformation block 455 which yields flight path 
rates and eventually body rates. Body angles Y el and and the velocity and 
acceleration vectors of the missile are functions external to the control system. 
These are not part of the actual mechanism of the guidance system and only 
facilitate to derive a representative input to the internal reference unit (IRU) and 
the radome. Periodically adaptive all angle guidance and control system 460 is 
represented as blocks in Figure 1, vectored LOS reconstruction 461 which supplies 
I G t * which is reconstructed LOS vector measurement of the range between the 



-27- 

target and the missile relative to a guidance frame. Block 462 creates an estimate 
of the range R tm from the target which when multiplied by the reconstructed LOS 
vector becomes a range vector R 0 *^ in the guidance frame. This input is received 
by the periodically adaptive object state estimators 464 and all angle guidance and 
control block 465, both of which receive input from IRU located in the missile. 
Other such systems with radome, seekers and autopilots can be used in conjunction 
with the periodically adaptive guidance as described above. 
[0058] Vectored LOS reconstruction 461 with a gimbaled seeker as well as 
with an electrically steered seeker, the periodic adaptive control and periodically 
adaptive target state estimation are now discussed in greater detail: 

2. Vectored LOS Reconstruction 

[0059] A simplified diagram of a means for creating a vectored LOS in a 
guidance frame from an object position parameter is presented in Figure 3. In 
particular Figure 3 shows vectored LOS reconstruction with a gimbaled seeker. 
The inputs to the vectored LOS reconstruction can include the seeker track errors 
e* el and e*^, seeker gyro rates co* el and co*^ and the IRU missile body to inertial 
navigational frame transformation matrix C n b , as well as the seeker Gimbal angles 
G* pitch and G* yaw . From the G* pitch and G* yaw , a transformation vector C B S from the 
missile body frame to the seeker bore sight frame is constructed in block 501. 
This transformation matrix is applied in block 502 along with the missile body to 
inertial navigation frame transformation matrix C G B and an identity matrix I which 
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equals the transformation matrix of the seeker to the inertial guidance frame at the 
initial condition of t=0. The output of block 502 is a transformation matrix at 
t=0 of the inertial navigation frame to the inertial guidance frame. 
[0060] The transformation matrix C g n (0) operates on the C n b matrix in block 
503 resulting in a transformation matrix from the missile body to the inertial 
guidance frame C G B . This transformation matrix operates on the transformation 
matrix C B S in block 504 to form matrix C G S which is the seeker bore sight frame to 
the guidance frame transformation matrix. The matrix C G S operates on the seeker 
gyro rates co* el and co*^ in block 505 forming the inertial seeker boresight direction 
cosine rate vector I G * S . The matrix C G S also operates on the seeker track errors e* e i 
and z\ forming a target end point to seeker boresight error vector e°V 
[0061] The direction cosine rate factor I e * s undergoes inertially referenced 
trapezoidal integration in block 520. The integration within blocks 506 and 512 
are accomplished using Z transforms in that I G * S is not a continuous function and 
thus lends itself more readily to the Z transformations. The output of the 
trapezoidal integration is the inertial seeker boresight vector I G * S which is 
combined with the error vector e 0 ^ at combiner 514. The resultant vector is 
normalized in block 515 to form a target aim point LOS vector I G * t with the 
magnitude of 1 , in other terms a unit vector. Referring to Figure 1 at block 462, a 
range estimate independently calculated is then applied to the LOS vector I G * t to 
obtain the range vector R 0 *^. The range vector and the IRU of the interceptor are 
used in the periodically adapted estimation and control. 
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[0062] Alternative means for creating a vectored LOS in a guidance frame from 
an object position parameter using vectored LOS reconstruction is shown in Figure 
4 as a simplified schematic wherein the seeker is electronically steered. The 
vectored LOS reconstruction with an electronically steered seeker is similar to the 
LOS reconstruction with the gimbal seeker, thus only differences of note will be 
discussed. Whereas in the gimbal seeker, gyro rates and gimbal angles are 
supplied as inputs, with an electronically steered seeker, these inputs are generated 
within the LOS reconstruction. However, the uncompensated seeker track errors 
remain as inputs to the electronically steered seeker LOS reconstruction. 
[0063] The seeker beam steering command space rates are generated in a 
control loop in which these uncompensated seeker track errors e el and are 
combined with the transformation matrix C g sc between the seeker beam steering 
command and the guidance frame in block 602. This track compensation with rate 
aiding results in the steering command space rates co e i and cv I G SC is integrated in 
an inertially referenced trapezoid integration process 610 resulting in I G SC the 
inertial seeker beam's steering command vector. The vector I° 5C is then added at 
612 to the target aim point to seeker beam center error vector e 0 ^, normalized and 
results in the target aim point LOS vector I G * t . The seeker beam steering 
command vector I G SC is further added at 614 to / ° sc operated on by [T s + T delay ] in 
box 620. The resultant vector is normalized 615, integrated 616, and multiplied 
617 by the transformation vector C G b which results in the seeker beam steering 
command vector in missile body frame I B fc which functions on the steering 
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command gimbal angles in block 618 generating the pitch and yaw angles G pitch 
and G yaw . Other functions shown in Figure 4 perform similar to their counterparts 
in the gimbaled seeker and are references as such. 

3. Periodically Adaptive Guidance 

[0064] Means for generating a guidance command signal using estimated object 
states utilizing an exemplary periodically adaptive guidance and control system is 
shown on Figure 5. The inputs to the guidance system are the reconstructed LOS 
vector I G * t , the seeker range and range rate measurements are R* ta , R * m , and IRU 
missile acceleration vector in an inertial navigation frame A n * m . In block 701, 
from the range and range rate measurements an estimate is made of range and 
range rate. The range rate k Tm is operated on by block 702 resulting in V c , the 
estimated closing velocity. The estimated range, R Tm is divided by V c in block 
703 to obtain the estimated time-to-go, or time to intercept t go . The means for 
generating a guidance command signal uses sequential line-of-sight (LOS) vectors 
in the guidance frame, in that the estimated range R Tm is multiplied in block 704 
by the LOS unit vector to obtain the range vector R 0 *^. The IRU missile 
acceleration vector in inertial navigation frame, A N * m , is applied to a 
transformation matrix in block 705 in order to obtain the missile acceleration 
vector A G * M in the inertial guidance frame. The periodically adaptive target state 
estimation in block 800 receives R^ and A G * m as inputs. As a result of the 
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periodically adaptive target state estimator 800, the means for generating a 
guidance command signal is periodically adaptive. 

[0065] The target state estimation result in four target state vectors as well as a 
m Ul model probability weighing vector W G m . The target state variables can be the 
estimated target to missile range vector Rj m , the estimated target to missile 
velocity vector Vj m and two target acceleration state vectors Aj m and Aj m . The 
weighing vector is represented as W G m . The estimated range vector is divided by 
the estimated range in block 706 to obtain a estimated LOS vector Ij , the 
estimated LOS vector is reduced into its constituents parts forming an estimated 
J G X vector in the guidance frame and an estimated K°_ T vector, both vectors being 
mutually perpendicular and formed in blocks 707 and 708 respectively. The 
estimated target to missile velocity vector is crossed by the estimated target to 
missile range vector and divided by the square of the scalar estimated range in 
block 709 thus resulting in an estimated target line-of-sight angular rate vector 
I g which is multiplied by the estimated closing velocity scalar V c in block 710. 
The target acceleration state vectors Aj m and Aj are multiplied by the m Ul model 
periodically adaptive control vector gains C G 3ra and C G 3 . m , and summed over the 
number of banks in block 720, the control vector gains are determined from the 
m* model probability Weight vector W G M in block 71 1 and 712. 0 13m (t g0 ) and 
$i4m(t g o) are derived in the previously mentioned master's thesis. 
[0066] The gains and navigation ratios associated with the exemplary 
periodically Adaptive guidance control system are shown in Table 2. The gains 
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C G 3m and C G 3 . m are used as the C 3 gains disclosed in prior art systems discussed 
previously. Laws C PA and D pa are associated with laws C and D respectively, but 
contained periodically adaptive gains. 
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t = Estimated-Time-to-Go [sec] 



^ 13m^ ^ 14m^J = Harmonican y Balanced Multiple Model PCA 
State Transition Matrix Elements 



-34- 

[0067] The periodically Adaptive Guidance Control system can include, as seen 
in exemplary embodiments, a gain C 4 in block 722 to compensate for autopilot 
lag, but is not required. 

[0068] The result of the periodically adaptive control calculation of block 720 
is an estimated zero effort miss vector ZEM f m divided by the time to go squared 
(i go ) 2 . The missile acceleration vector in an inertial guidance frame A G * m is 
subject to a gain in block 722 of C 4 , similar to classical optimal control theories 
discussed previously- The LOS component vectors J G t , K G , the quantity V c and 
j G , the Zero effort miss vector ZEMtm divided by (^ 0 ) 2 and the missile 
acceleration vector multiplied by the gain C 4 combine in block 725 to obtain a 
missile guidance acceleration command vector N G C in the inertial guidance frame 
which is perpendicular to the LOS. A close look at the equation terms in block 
725, reveals the means for generating a guidance command signal comprises a 
proportional navigation controller, a classical optimal controller and an augmented 
proportional navigational controller. The coefficient V C 'X ° originates from 
proportional navigational control theory. The coefficient C 4 A G * m represents the 
acceleration command that compensates for autopilot lag which is synonymous 
with classical Optimal Proportional Navigation Guidance. The ZEM tm vector 
divided by time-to-go squared (i go ) 2 is similar to the augmented acceleration 
command vector and is the predominate periodic adaptive guidance portion. The 
resultant missile guidance acceleration command signal is a vector that is 
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perpendicular to the LOS vector / f thus the significance of the "skid to turn" 
missile maneuver term can be seen. 

4. Target State Estimation 

[0069] Means for estimating object states according to an exemplary 
periodically adaptive target state estimation system referred to here before in 
general terms is shown in greater detail in Figure 6. The inputs to the periodically 
adaptive target state estimation system are the reconstructed LOS vector, R G * Un , a 
seeker measurement error of variance matrix transform from seeker boresight 
frame to inertial guidance frame [a 2 r ] G and the IRU missile acceleration vector n 
inertial guidance frame A G * m . The target state estimation is carried out in each of 
the three axis independently. To be periodically adaptive, the means for generating 
a guidance command signal has at least one Kalman filter bank, to be periodically 
adaptive, the bank is harmonically balanced, and each of the at least one Kalman 
filter banks 810a, 810b and/or 810c is associated with a respective axis x,y and/or 
z in the guidance frame G. Although Figure 6 shows the estimation of these target 
state variables for each of the three axis, for description purposes only the x axis 
will be discussed as the process is identical for each axis in the guidance frame. 
[0070] The harmonically balanced Kalman filter bank 810a receives the X 
component of R 0 *^ and designated as R GX * , the expected X axis target 
acceleration variance o\ x and the X component of the seeker measurement error 
variance matrix in the inertial guidance frame o 2 RX and the IRU missile 
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acceleration vector and inertial frame A Gx m . In generating a guidance command 
signal, the Harmonically Balanced Kalman filter Bank 810 is a means for 
generating a set of probability weights. An output of the harmonically balanced 
Kalman filter bank 810 is the m* model X axis probability weights W Gx m and the 
m Ul model X axis state estimates X m which comprises the vector/?^ , V^ m » 
and i GX . The means for generating a guidance command signal creates an 
estimated object to device range vector, an estimated object to device velocity 
vector, an estimated object acceleration vector and an estimated object acceleration 
rate vector, these vectors are representative of the object or target states shown in 
Figure 5. The summation in block 81 1 of the product of W Gx in and X m where m 
= 1 to number of banks, results in the X axis state estimate X which includes 
R™* , V™ , A r X , and • gx. The combination of all three dimensions, results in 

TM 7 1M ' 1 7 A. j 

an output of the periodically adaptive estimated states of an estimated 
target- to-missile range vector R G TM , an estimated target-to-missile velocity vector 
Vj M and the m Ul model estimated target acceleration state vectors Aj and ^ g also 
containing three components associated with the axis of the guidance frame. 
Additionally, the m Ul model probability weighing vector W G m is provided. It is 
noteworthy to mention that each axis state estimates can be generated 
independently and congruently. 



-37- 

5. Harmonically Balanced Kalman Filter Banks 

[0071] Kalman filtering theory models the system using Gaussian distributed 
random noise disturbances. The Kalman filter is a real-time computer algorithm 
that estimates system states with a recursive two step process. Kalman filters 
predict system states expected at the next measurement input time using a priori 
knowledge of the missile - target dynamic model, and it corrects the predicted 
object states with noise corrupted measurements of the system. 
[0072] A simplified block diagram of an exemplary harmonically balanced 
Kalman filter in the X axis is presented in Figure 7. Other harmonically balanced 
Kalman filter banks for the Y and Z axis are performed similarly. 
[0073] A Bayesian hypothesis probability calculator is represented in the 
diagram as area 910. The inputs to the Bayesian hypothesis calculator are 
o resX m (t K ) which is the multiple model X axis residual variances which is matrix 
of constants, the derivation of which is well known to those skilled in the art. Res 
X m (t K + ), the multiple X axis residuals also serves as inputs. An m dl model 
Kalman filter with the "correct" system dynamics hypothesis model can be 
expected to have minimized residuals compared to other "mismatched model" m m 
model Kalman filters. From these two inputs, a multiple model of Gaussian 
probability density distribution is created at block 917. The measurement 
probability density of the m* filter in the multiple model bank is a Gaussian 
distribution, in that the Kalman filter theory depends on the system modeling with 
Gaussian probability density distributions. 
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[0074] Throughout Figure 7, (t K ~) is the unconditioned result of a prediction 
estimate while (t K + ) is the conditioned prediction estimate. The multiple model 
Gaussian probability density 917 is multiplied by the multiple model X axis 
probability weights in block 911. The result is then normalized in block 912 by 
dividing the sum of the product of the m* model Gaussian probability density with 
the m Ul model X axis probability weight from one to the number of filter banks 
(NMB). The result is a multiple model X axis probability weight that has been 
conditioned. 

[0075] The multiple model W Xm probability weights are integrated in block 
913 using a Z transformation with W Xm (0) being the initial condition for the 
transformation. The m m model X axis probability at (t K _ t ) which is unconditioned 
is multiplied by the transition probability matrix T MEQMP in block 914. T MEQMP is 
the transition probability of the true system dynamical parameters not changing 
during sample time interval and is set at 99.9%. T MNEMP is the transition 
probability of true system dynamical parameters changing during a sample time 
interval and includes all values where T MNEMP is not equal to T MEQMP . The 
formulation of this matrix and it's use in Bayesian hypothesis probability is known 
and documented as the adaptive learning theory of Moose and Wang entitled "An 
adaptive estimator with learning for a plant containing Semi-Markov switching 
parameters, " by R.L. Moose and P.P. Wang, IEEE Transactions on Systems, 
Man, and Cybernetics, vol. 3, (May 1973): 277-281, which is herein incorporated 
by reference. The matrix is summed with T MNEMP in block 915 which results in 
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the multiple model X axis probability weight W Xm (t K ~) which is the unconditioned 
result of prediction estimates. The use of the transitional probability matrix allows 
the bank of harmonically balanced Kalman filters to move from weights =0 to 
weights > 0. 

[0076] The X axis missile acceleration measurement U x (t K ) is multiplied by L 
[- T s 2 /2, -T s ,0 ,0] T in block 929 and added to the.PCA state transition matrix 
^ > m PCA ( T s) where T s is the sample time in block 930 to form the deterministic input 
XJt k +\) . The deterministic input is integrated to form XJt k -) through a Z 
transformation in block 931. The initial input into the Z transform is the initial 
condition of X at T=0. The result is multiplied by the H vector which is the 
system measurement matrix shown in block 932 and subtracted from the X axis 
range vector X* (t k ) in block 901 . The result of this addition are multiple model X 
axis residuals, res Xm (t K ), which continues the loop as input to the multiple model 
Kalman probability function of block 917 and are summed in block 934 along with 
a multiple model X axis Kalman gain matrix K Xin (t K ). This combination added 
to X^t K -) at 933 supplies X m (t k + ) which is a multiple model X axis state 
estimate. 

[0077] The expected X axis target acceleration variance o 2 ax is multiplied by 2 
A Tm which reflects the plant noise and r m r m TPCA (T s ) in block 921 and added to the 
covariance propagator in block 922. The result is the multiple model X axis 
covariance matrices P Xm (t K ") which is thereafter multiplied by the transpose of the 
transition matrix H and divided by the residual variances o^ 2 ^ (t K ") in block 921. 
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1=1: 



The result is the multiple model X axis Kalman gain matrix K Xm (t^ which is 
applied to block 925 which conditions the multiple model X axis covariance 
matrices. The conditioned multiple model X axis covariance matrices are 
submitted to block 926 to limit the cross variance, so that the correlation of 
coefficients is < 1. Covariance limiting is required to eliminate singularities in 
the multiple model X axis covariance matrices. The result is integrated in block 
927, by a Z transformation, and multiplied in block 922 to form covariance 
propagation. 

[0078] When an initial condition is true as in, block 935, no residue is created 
as the path is broken, however the propagation of the multiple model X axis 
covariance matrices P Xm continues for multiple time periods T s . The result of this 
process is X m {t K +) which is a multiple model X axis state estimate. The 
harmonically balanced multiple model PCA noise covariance input matrix r m r m T 
«* (T s ) is digitally precomputed from the equation (24). 



[0079] It is important to note that T s is a fixed sample time in seconds. Thus 
the Bayesian hypothesis estimator the covariance propagator, the Kalman matrices 
and the multiple model target state estimates result in the multiple model X axis 
state estimates X m (t k ) , which is conditioned. 




(24) 



o 



where G m = [0 0 1 (^ Nm ~2X Tm )f 
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[0080] The expected X-axis target acceleration variance a ax 2 is shown in the 
top, left hand corner of Figure 7. This is a constant that is selected using the 
generalized knowledge of the target and is constant within the velocity frame. 
Multiple a ax 2 can be used depending upon the type and nature of the target 
encountered. An independent system or operator input can be used to evaluate and 
select from a plurality of established a ax 2 based on a priori knowledge of the 
specific targets. 

[0081] Figure 8a is a representation of a flight path 866a of a missile 860a 
intercepting a target 850a, executing a barrel roll maneuver flight path 855a, using 
an exemplary embodiment of vectored LOS reconstruction and Classical Optimal 
Guidance Control. The figure demonstrates the guidance control systems ability to 
engage in high boresight angle engagements. The Classical Optimal Guidance 
Control results in the missile mimicking the targets flight path with an expenditure 
of additional energy. The Figure also contains an insert 870a which is the target 
acceleration power spectral density plot. 

[0082] The coverage of classical optimal control law dynamics is shown where 
the object (target) bandwidth X T is one radian per second and the normalized target 
acceleration has a frequency of o> T = 0, both of which are assumed from a priori 
knowledge of the target. The barrel roll maneuver at co T = 3, is outside the 
response capability of the classical optimal control law dynamic, the consequence, 
as shown in Figure 8a, being the intercepting device 860a mimics the target's path 
855a as seen the Figure. . 
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[0083] Figure 8b represents a flight path 866b of a missile 860b intercepting a 
target 850b, executing a similar maneuver along flight path 855b, using an 
exemplary embodiment of Periodically Adaptive Guidance Control. The figure 
demonstrates the Periodically adaptive guidance control system's ability to adapt to 
the maneuver frequency of the target and minimize the necessity, to match the 
periodic maneuvers of the target. The periodically adaptive guidance control 
system can reduce energy consumption during intercept as well as required missile 
agility. The Figure also contains an insert 870b which is the target acceleration 
power spectral density plot. 

[0084] The periodically adaptive control law dynamics as seen in Figure 870b, 
encompasses the barrel roll maneuver frequency o> T = 3 within a small bandwidth 
.25 rads per second. The periodically adaptive capability allows the system to 
adapt to varying target acceleration frequencies. Thus the control system can 
respond robustly to the periodic barrel roll or other periodic maneuvers. The its 
resultant missile path shown as 866b. 

[0085] Figure 9a shows a periodically adaptive control law dynamic of an 
exemplary embodiment. This Figure is reproduced from Figure 8b such that the 
adaptive dynamic nature of the present invention can be illustrated when view 
alongside Figure 9b- Figure 9d. 

[0086] Figure 9b is a representation of the bandwidth response capability for a 
constant maneuver target where the maneuver frequency co T = 0. The classical 
optimal control dynamics can encompasses the constant acceleration maneuver 
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dynamics of the target, since it is within the assumed bandwidth. Figure 9b shows 
that the periodically adaptive control dynamics can effectively encompass the 
target acceleration frequency of the constant maneuver target. 
[0087] An overlap in the classic optimal control law dynamics response 
bandwidth can occur for a jinking target. However, as the bandwidth of the target 
maneuver jinking dynamics increases, the effectiveness of the classical optimal 
control dynamics is reduced. Figure 9c is a target acceleration power spectral 
density plot for periodically adaptive control against a time varying jinking target. 
Figure 9c demonstrates that the periodically adaptive control dynamics response 
bandwidth adapts to encompass the bandwidth and frequency of the target. The 
same is true of Figure 9d where a constant velocity target is presented. 
[0088] As is evident in Figures 9a-9d, the periodically adaptive control 
dynamics response bandwidth of the present invention readily adapts to target 
maneuvers that are periodic with a frequency greater than zero. In the case of 
constant acceleration, constant maneuver and jinking maneuver, the periodically 
adaptive control system demonstrates a capability at least equal to those of a 
classical optimal guidance system. Therefore, in addition to handling target 
maneuvers in which the classical optimal guidance system is well formed, 
exemplary embodiments of the present invention can handle a periodic maneuver, 
such as the barrel roll maneuver used by aircraft pilots for evasive maneuvers, or 
the tumbling of intercontinental ballistic missiles upon their reentry. 
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[0089] An exemplary embodiment utilizing alternate means for determining a 
guidance command using estimated object state with classical optimal guidance is 
presented in Figure 10. Classical optimal all angle guidance shares similar 
components and operations with exemplary embodiments of periodically adaptive 
all angle guidance system discussed with respect to Figure 5 and are referenced 
accordingly, with the notable exception of periodically adaptive object state 
estimation and periodically adaptive control. 

[0090] The inputs to the guidance system are the reconstructed LOS vector I G * t , 
the seeker range and range rate measurements are R* to , R* m , IRU missile 
acceleration vector in an inertial navigation frame A n * m . In block 701, from the 
range and range rate measurements an estimate is made of range and range rate. 
The range rate is operated on by block 702 resulting in V c , the estimated 
closing velocity. The estimated range, R m is divided by V c in block 703 to obtain 
the estimated time-to-go, or time to intercept t go . The estimated range R tm is 
multiplied in block 704 by the LOS unit vector to obtain the range vector R 0 *^. 
The IRU missile acceleration vector in inertial navigation frame A N * m is applied to 
the transformation matrix in block 1405 in order to obtain the missile acceleration 
vector A G * M in the inertial guidance frame. The classical optimal target state 
estimation function in block 1400 is well known in the art and it particulars are not 
further developed. The classical optimal target state estimation function 1400 
receives R 8 *^ and A G * m as inputs. 
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[0091] In an exemplary embodiment, means for estimating object states result 

in three target state vectors. The target state variables are the estimated target to 

missile range vector R f m , the estimated target to missile velocity vector V f m and 

the estimated target acceleration state vector The estimated range vector and 

velocity vectors are processed in the same manners previously discussed regarding 

Figure 5, resulting in an estimated J G vector in the guidance frame J f and an 

estimated K°. vector K f , and an estimated target LOS angular rate vector 

• G which is multiplied by the estimated closing velocity scalar V c in block 710. 
A 

The target acceleration state vector A t ° is multiplied by gain G> at block 1420 
where Q is a function of A x and t g0 ,as shown in equation (25) where X T is a 
predetermined target bandwidth based on priori knowledge. 



e~ Xj + X r -l 

C 3 = (25) 



where Xj " XjX ' 8 ° 

X T = Target Maneuver Bandwidth 

[0092] The result of this gain in block 1420 is an estimated zero effort miss 
vector ZEM^ divided by the time-to-go squared i^ Q . The missile acceleration 
vector in an inertial guidance frame A G * m is subject to a classical optimal control 
theory gain of C 4 in block 722 as discussed previously. The LOS component 
vectors J f , K f, the quantity V c and the Zero effort miss control vector 
ZEM %m divided by i* 0 and the missile acceleration vector multiplied by the gain 
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C 4 combine in block 725 to obtain a missile guidance acceleration command vector 
N G C in the inertial guidance frame which is perpendicular to the LOS. 
[0093] Figure 1 1 shows an exemplary embodiment of the present invention 
applied in a fire control system. The system containing means for generating a 
signal representing a predicted position of the object from object position 
parameters relative to a guidance frame and a periodically adaptive estimated 
object state produced in the guidance frame using the object position parameters; 
and, means for transmitting the signal to an on-board guidance control of the 
device. 

[0094] The fire control platform 1500 tracks the target 1520 through a radar 
system 1501 in addition also tracks the interceptor missile 1510. The fire control 
platform can be a fixed ground station, mobile station an aircraft, ship, or space 
based platform. The missile fire control platform includes a radar processing unit 
1503 that determines range and directional information of the target and of the 
interceptor missile. These values are fed into periodically adaptive object state 
estimators 1504 which create a predicted intercept position (PIP). The PIP is then 
transmitted through a transmitter 1505 via an antenna 1502 to a receiver on the 
interceptor missile 1511 where the missile generates an acceleration command 
based upon the PIP. Whereas previously exemplary embodiments of the invention 
are successfully used for terminal guidance, missile file control guidance is 
particularly suited toward mid-course guidance. 
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[0095] Interceptor missiles due to constraints of energy and weight have a 
narrow acquisition band which is well suited to close in targets, however a missile 
fire control radar system, where power supply and size of the radar system is not a 
limiting factor, can more readily acquire the targets. The missile fire control 
system can direct the interceptor missile until the interceptor missile can acquire 
the target with its own tracking system, or direct the missile through termination. 
[0096] Figure 12 shows an exemplary embodiment of periodically adaptive 
missile fire control system. A missile fire control platform 1600 obtains the range 
and range rate measurements from it's tracking radar and forms range and range 
rate estimation in block 1601. These inputs can be the measured range of the 
target to the platform R* TP , the range rate of from the target to the platform R* p , 
the range of the missile to the platform R^ P and the range rate to the platform 
R^ P > The estimating closing velocity V,can be found by subtracting the range 
rate of the missile to platform from the range rate of the target to platform. From 
this the closing velocity is found. Also in block 1603 the range of the target to 
missile is divided by the closing velocity to obtain an estimated time-to-go t go . 
The range vector determined from the radar inputs is given as R E XP wherein E is in 
the earth frame. 

[0097] The periodically adaptive target state estimation creates an estimated 
target to platform range vector R^ p , an estimated target to platform velocity vector 
Vj P and multiple model target acceleration state vectors Aj m and Aj m as well as 
the multiple model probability weighting vectors W ra In the periodic control box 
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of 1605 the multiple model estimated target acceleration state vectors and the 
multiple model probability weighing vectors are multiplied and summed over the 
number of Kalman banks in each respective axis. The result is a zero effort miss 
vector in the earth frame. The output is a predictive intercept position (PIP) 
formed in block 1606. The PIP is simply the range Rj P plus the velocity vf P 
multiplied by the time-to-go i plus the zero effort miss vector ZEM^ along 
with the estimated time-to-go. The PIP is an prediction of the future position of 
the object at (T + 1 ), where T is real time. . 

[0098] This information is transmitted from the missile fire control transmitter 
1607. The method of transmission to the missile interceptor 1610 can include RF, 
wireless or wired means. The uplink is received in receiver 1611 which is fed into 
the midcourse guidance system 1612 which uses the uplink signal and information 
form an inertial navigation system 1613 to generate a solution and a missle 
guidance command enabling the missile to intercept the target. The missile 
guidance command is relayed to the autopilot 1614 and ultimately to the thruster 
or other control apparatus. 

[0099] Because of the relatively non-dynamic relationship between the missile 
fire control platform and the target small angle approximations can suffice in 
creating the PIP and reducing the necessity to create a reconstructed vectored 
LOS. 

[00100] Although the embodiments are directed to an air-to-air missiles, the 
application of the claimed subject matter to air-to-surface, surface-to-surface, and 



